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Abstract 

We study ergodic properties of a family of traffic maps acting in the space of 
bi-infinite sequences of real numbers. The corresponding dynamics mimics the mo- 
tion of vehicles in a simple traffic flow, which explains the name. Using connections 
to topological Markov chains we obtain nontrivial invariant measures, prove their 
stochastic stability, and calculate the topological entropy. Technically these results 
in the deterministic setting are related to the construction of measures of maximal 
entropy via measures uniformly distributed on periodic points of a given period, 
while in the random setting we directly construct (spatially) Markov invariant mea- 
sures. In distinction to conventional results the limiting measures in non-lattice 
case are non-ergodic. Average velocity of individual "vehicles" as a function of their 
density and its stochastic stability is studied as well. 
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1 Introduction 



Speaking about stochastic stability of a dynamical system one means (see e.g. [24, 4]) 
that the most important statistical quantities related to the dynamics (e.g. Sinai-Bowen- 
Ruelle invariant measures) depend continuously on the addition of a small amount of true 
random noise. We study ergodic properties of a family of traffic maps acting in the space 
of bi-infinite sequences. The corresponding dynamics mimics the motion of vehicles in a 
simple traffic flow, which explains the name. Our aim is to construct nontrivial invariant 
measures of the traffic maps and to show that they are stable with respect to "natural" 
random perturbations - by which we mean that the motion of individual vehicles are 
performed with a probability p. Then the case p = 1 corresponds to the deterministic 
traffic map, while situations with p < 1 may be considered as random perturbations. 
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We introduce the following notation. Under an admissible configuration at time 
t e Z_|_ U {0} we mean an ordered countable set of particles (balls) of radius r > 0, centers 
of which are located at the points := {xl )iez C R C such that 

xl + r< - r. 

The set of all admissible configurations we denote hy X — X{r,'R). By v > denote the 
maximal possible movement of a single particle per unit time, i.e. 

< xl+^ -xl<v. 

The parameter p G (0, 1] stands for the probability of movement for individual particles. 
A single particle performs a totally asymmetric random walk (which explains why the 
processes of this sort are often called TASEP - Totally Asymmetric Exclusion Process) 
with jumps of size v, occurring with probability p, until its motion does not interfere with 
the motions of other particles. 

Exclusion processes (collective random walks of countable collections of particles with 
hard core interactions) introduced by Frank Spitzer in 1970 appear naturally in a broad 
list of scientific applications starting from various models of traffic flows [20, 17, 12, 
5, 6], molecular motors and protein synthesis in biology, surface growth or percolation 
processes in physics (see [21, 9] for a review), and up to the analysis of Young diagrams in 
Representation Theory [11]. Continuous time versions of these processes are reasonably 
well understood (see e.g. [19] for a general account and [1, 2, 13, 15] for recent results). 
The main difficulty in the analysis of discrete time versions of exclusion processes is that 
an arbitrary (infinite) number of particle interactions may happen simultaneously. To 
overcome this difficulty one needs to develop principally new approaches. Note that in 
the one-dimensional setting under consideration the basic restriction is the preservation 
of the order of particles, i.e. the particles can not overtake each other under dynamics. 

We consider discrete time Markov processes 7r(p, i',r, R) acting in the space of con- 
figurations X — X{r, R) and the dynamics of individual particles in the configuration is 
defined by the relation 

t+i _ (mm{xl + V, xj^i — 2r} with probability p 

' \xl with probability 1 — p ' \ ■ ■> 
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Figure 1: Exclusion process in continuum 

Consider three, at first glance, very difi'erent types of exclusion processes satisfying 
the relation (1.1). Processes of type 1 act on a lattice R = Z^,r = 1/2, ^; e One 
lattice site may be occupied by at most one particle. Models of this type are widely used 
to describe the motion of vehicles on a single-lane road (see for example [20, 22]). 

Processes of type 2 also act on a lattice R = Z^, v G Tj^, but r = 0. The fundamental 
difference is that a lattice site can be occupied by an arbitrary number of particles. Mod- 
els of this sort with a continuous-time are called zero-range processes (see, for example. 



[14]) and it is convenient to use them to simulate a communication line, in which parti- 
cles represent equal packets of information, waiting in queues to communication servers, 
located at sites of the lattice R. In terms of quantum statistical mechanics processes of 
type 1 and 2 are related as interacting Fermi gas and free Bose gas. 

Processes of type 3, which are a special case of exclusion type processes introduced in 
[7], act not on a lattice but on the continuous space R = IR^,r>0,i'e R^. It will be 
shown (Lemma 1) that for r = 1/2 these processes contain all realizations of the processes 
of the 1st type, and for r = all realizations of the processes the 2nd type. This allows one 
to obtain simultaneously analytical results for all cases. The consideration of exclusion 
type processes in continuum not only simplifies the analysis but due to the presence of 
additional symmetries (absent in lattice versions) it offers a possibility to obtain new 
results about the lattice cases unavailable otherwise. 

A density of a configuration (the number of particles per unit length) is defined as 

q(x^) := lim n/{xl^_^-xl), 

n— >-oo 

if the latter limit makes sense. ^ As we shall see the density is the first integral of the 
processes under study and thus the phase space is foliated naturally by invariant subsets 
X^^^ consisting of admissible configurations of density g. 

The formula (1.1) describes the processes through dynamics of configurations of or- 
dered particles. This is convenient for the analysis of the particles motion, but does not 
allow the study of invariant measures (stationary distributions). To this end we consider 
a modification of the process under study in which the particles are indistinguishable from 
each other. Speaking about invariant measures we always refer to the invariant measures 
of this modification. The set of probability n{p,v, r, R)-invariant measures we denote by 

•M.p,v,r,Ti- 

To date, the mathematical description of invariant measures for the deterministic pro- 
cesses under study is essentially absent, and the only classification result [3] gives just 
a formal description of invariant measures for the simplest lattice process 7r(p = 1,^ = 
l,r e {0, 1/2}, Z) and says nothing even about the existence of nontrivial (non-atomic) 
invariant measures. Note that the processes Tr{p — l,v,r,Ii) possess periodic trajectories 
of all possible periods and hence atomic measures supported by these trajectories. In 
Theorems 1,4,5 for each density g we prove the existence of a non-atomic invariant mea- 
sure iig supported by configurations of the given density. Moreover we give an explicit 
construction in terms of (spatially) Markov measures related so subshifts of finite type. 

Having a large number of invariant measures it is important to distinguish "physically 
relevant" ones. In the case of low dimensional dynamical systems one often uses for 
this purpose the concept of Sinai-Bowen-Ruelle (SBR) measures. Roughly speaking the 
latter means that for a reasonably large family of "good" initial probabilistic measures 
(say all absolutely continuous measures) their images under dynamics converge in Cesaro 
means to the same SBR measure. In the present infinite-dimensional setting the choice 
of "good" initial measures is not obvious and the control over the convergence of their 
images is not available at present. An alternative approach consists in the analysis of 
stochastic stability of invariant measures. In a number of cases it has been shown (see, 
e.g., [4]) that SBR measures are exactly those that are stochastically stable. 

^The one-sidedness of the growth of segments {xq,xIi__i) is due to the fact that aU particles move in 
the same direction. The definition of the density in the general case is more complicated (see e.g. [7]). 



As we already noted a natural choice of random perturbations for the deterministic 
processes 7r(p = l,v,r,H) are random processes 7r(p < l,v,r,Ii). Here again in the 
discrete time case only very partial results about invariant measures are known in the 
literature. Nevertheless the simplest lattice case 7r(p < 1, v = 1, r = 1/2, Z) and especially 
its much simpler continuous time version was intensively studied in physics literature from 
this point of view (see [1, 13, 14, 15, 16, 20, 22, 23]). In particular, in [22, 23] using 
a mean field approximation and an interesting and nontrivial combinatorial argument 
(apparently not quite complete without the exact analysis of the limit construction) the 
authors have found necessary conditions (equivalent to our formula (4.2) for the existence 
of nontrivial invariant measures. On the other hand, the next step - the proof that the 
constructed measures are indeed invariant under dynamics was not explicitly addressed. 
The construction proposed in the proof of Lemma 4 provides not only a complete proof of 
this statement but also in its first few lines gives the crucial formula (4.2) without any limit 
transitions. By means of a completely different approach M.Kanai [18] has studied the 
same process on a finite discrete ring (instead of the infinite lattice) and obtained an exact 
rather complicated formula for the invariant measure as a function of the ring's length 
L. When L goes to infinity this invariant measure converges to the limit derived in [22], 
which is of independent interest despite a number of restrictions assumed in this paper. 
The approaches elaborated in [22, 23, 18] use in one way or another properties of invariant 
measures of the process 7r{p,v = l,r = 1/2,{1,2,...,L}) on a finite ring. Actually in 
[22, 23] the Markov structure of the invariant measure was discovered (without the proof 
of its existence) while in [18] the unique invariant measure is no longer Markovian and 
becomes Markovian only in the limit as L — > cxo. A coupling construction proposed by 
L. Gray [16] allows one to show that for a fixed density g the Markov invariant measure 
is ergodic and unique among translationally invariant measures. 

It is worth noting that the above mentioned lattice constructions cannot be extended 
to the case of long jumps {v > 1). To overcome this difficulty we first study the continuous 
setting to obtain the invariant measure for the process 7r{p < l,v = l,r = 1/2, R) and 
then extend this result for all values of r using a wide set of symmetries available in R 
(but not in Z). After that we are able to return back to the lattice cases. 

Despite the fact that the constructions of invariant measures in deterministic and 
stochastic settings are rather different and neither of them can be applied in another 
setting, the proof of stochastic stability is relatively simple and is based on the explicit 
description of the invariant measures. 

Our main result about invariant measures is as follows. 

Theorem 1 The process 7r(p, f,r, R) for each < p < l,v > 0,r > and each particle 
density g G (0, l/(2r)] possesses a nontrivial invariant measure fjp^ supported by the set of 

configurations of density g and /^^ in the weak sense. 

One of the most important characteristics of the processes under consideration is an 
average velocity of particles during time t > 0: 

Vix,i,t) := {xl-x^,)/t. 

In [7] it was shown that under very general assumptions (including, in particular, dynamics 
of particles moving in opposite directions), which definitely hold true for all processes 
considered in the present work, the a.s. limit of the statistics V{x, i,t) as t ^ oo (if it 



exists) depends only on the density q of the configuration x (see Theorem 3) . Therefore, 
to calculate the average velocity 

V{x) = V{g{x)) := lim V{x,i,t) 

(where the convergence is considered in the almost sure sense) for a configuration of a given 
density it is enough to determine the latter for a specially selected initial configuration 
of the same density (for which the calculation is convenient). Note that the existence of 
the limit under question, not to mention the explicit formulas for it, previously has been 
proven only in the deterministic setting (i.e. for p — 1), see [7]. The following result 
gives explicit relations between the average velocity and other parameters of the process 
7r(p, V, r, M). 

Theorem 2 For the process 7r{p,v,r,R) \/v,q e R\, < r < 1/(2^), p e (0,1) with 
any admissible initial configuration of density g the average velocity V{g, r) is well defined 
and is calculated as follows: 

V{g,r) = {l-2rg) F(^-A_,0), 
1 + v/(TT^;^P^^4^] /(2^) ^min{l/^,^;}. (1.3) 
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Figure 2: Fundamental diagrams (average velocity against density) for the process 
7r(p, V, r = 0, M). 

In physics literature results similar to Theorem 2 are known for the simplest lattice 
TASEP process ^{p < l,v = l,r = 1/2, Z) (see [22, 23, 18]). Here again one needs to 
make comments about missing arguments of the same sort as in the analysis of invariant 
measures plus that the result should hold not only for "typical" initial configurations (see 
discussion below and in Section 5). 

The dynamical coupling construction developed in [7, 8] allows one to get complete 
information about the properties of average velocities in the deterministic setting (i.e. for 
p = 1). This construction does not require the study of (numerous) invariant measures of 
the process, but gives only conditional (upon their existence) results in the stochastic case 
(albeit under much broader assumptions about the process: the local velocities v = Vi 
are iid random variables). In the stochastic case one cannot avoid the analysis of invari- 
ant measures. Explicit expression for the invariant measure fi allows one to derive (by 
Birkhoff's ergodic theorem) the formula for the average velocity for /x-a.a. initial config- 
urations. To extend this result to a much more broad setting (formulated in Thereom 2) 



for all particle configurations for which densities are well defined, one needs to use met- 
ric properties of the process obtained earlier. Namely in [7] it has been shown that for 
all particle configurations of a given density the average velocity is the same. Therefore 
it is enough to calculate this statistics for a single suitable initial configuration (e.g. a 
configuration typical with respect to the measure /x). 

The paper is organized as follows. In Section 2 we discuss connections between the 
processes under study and briefly review known results related to particle densities and 
average velocities. Ergodic properties of the deterministic traffic map will be discussed 
in Section 3, which we shall finish by the calculation of the topological entropy of the 
traffic map in continuum. The true random setting (i.e. < p < 1) will be analyzed in 
Section 4. In Section 5 we shall study average particle velocities, obtain explicit relations 
between these statistics and parameters of the process, and prove their stochastic stability. 
Finally in Section 5.4 we shall discuss heterogeneous versions of the processes imder 
consideration (particles with different sizes in the same configuration and presence of 
randomly distributed static obstacles in space). 

2 Metric properties 

Here we shall discuss connections between the processes under study and briefly review 
known results related to particle densities and average velocities. 

Let Aj = Ai(x*) := xl_^_i — x\ — 2r stands for the distance between boundaries of the 
balls corresponding to subsequent particles in the conflguration a;*(see Fig. 1). We refer 
to Aj as the gap corresponding to the i-th particle in x*. If Aj(x*, z) <v we say that the 
i-th particle is blocked (meaning that its motion is blocked by the i + 1-th particle) at 
time t and free otherwise. By a cluster of particles in a conflguration G X we mean a 
locally maximal collection of consecutive blocked particles. 

To emphasize the dependence of various statistics on the parameter r > (ball's 
radius) we denote by x{r) and X(r, R) the conflgurations with balls of radius r and the 
corresponding space of admissible conflgurations. In the special case of point-particles 
(r = 0) the dependence on the radius will be omitted. 

We say that two particle processes whose dynamics is described by the relation (1.1) 
are statically coupled if all random choices related to the motion of particles with the same 
indices in these processes coincide. 

The following results show relations between the three types of exclusion processes 
(two lattice and one in continuum) introduced in the previous Section. 

Lemma 1 For any given v,p the following relations are valid: 

1. n{p,v e Z+,r = 1/2, Z) = n{p,v,r ^ 1/2, K) ifx^{r = 1/2) C Z; 

2. 7r(p, t; e Z+, r = 0, Z) = 7r(p, v,r^O, R) ifx^{r = 0) C Z; 

5. Vr > there exists an affine horneomorphism ip = '■ X(r, M) — )■ X(0,]R), such 
that 

</? o 7r(p, V, r, R) = 7r(p, t;, r = 0, R) o 

VM,f>0, Q < p < 1 a sub-lattice R„^^ := fZ + m is invariant with respect to the 
process 7r(p, t;, r = 0, R), i.e. x^ C ^ =^ C R^ ^ > 0. 



Proof. The only nontrivial statement here is the item (3). Observe that any two 
particle configurations x(r), x{r) having the same sequence of gaps A := {Aj} may be 
transformed to each other by a one-to-one map 

Xi(r) = ^r{xi{r)) := Xi(r) - 2i(r - r) Vi e Z. (2.1) 

This allows one to choose the homeomorphism ipr as follows 

{(Prix{r)))i := Xi{r) - 2ir Vi e Z. 

Now the conjugation between the processes having particles of different sizes follows from 
the equality between sequences of gaps A* between particles in the statically coupled 
processes 7r(p, v, r, R) and 7r(p, v, 0, R). □ 

The correspondence between densities and average particle velocities for statically 
coupled processes with configurations consisting of balls of different radiuses r ^ r > is 
summarized as follows. 

Lemma 2/7/ 

1. The density q{x*) is preserved by dynamics, i.e. q{x^) = g{x^^^) Vt > 0. 

2. Let configurations x(r) G X(r, M), r > and x{r) G X(r,M), r > have the same 
sequence of gaps {Aj}. Then 

Q{x{r)) - '^^^^^ 



l + 2(r-r)Q(x(r))' 

3. Let additionally the processes a;*(r) and i*(r) be statically coupled. Then \/i,t 

V{x{r),i,t)^V{x{r),i,t). 

Results obtained in [7] allow also to claim the coincidence of average velocities both 
for individual particles and for configurations having equal densities. 

Theorem 3 ([7]) For the process 7r(p, v, r, R) Vv, r G R, ^ > 0, p G (0, 1] the following 
claims are valid. Let x,y & X{r, H), g{x) — g{y) — g and let the average velocity V{y) be 
well defined. Then 

\V{x,z,t)-V{y,j,t)\'-^0 

for any i,j G Z. 

3 Ergodic properties: deterministic case (p=l) 

Despite the abundance of invariant measures supported by time-periodic trajectories (ex- 
isting for all periods) nothing was known previously about the existence of nontrivial 
invariant measures. Our main result in this direction in the case R := R is the following 
claim. 



Theorem 4 The deterministic process 7r(p = l,i',r, R) for each v > 0,r > and each 
particle density g e (0,1/ (2r)] possesses a non-atomic invariant measure /ig such that 

The proof of this result will be given step by step in Sections 3.1,3.2. In the simplest 
lattice case the construction of the invariant measure will be described explicitly in terms 
of the so called (spatially) Markov measures. The notion Markov measure here refers 
to the unique invariant measure for the Markov shift on {0, 1}^ with a non-degenerate 
transition matrix (pij), i,j G {0, 1}, Pij > 0, '^jPij = 1. In terms of configurations of 
particles, 1 corresponds to the presence of a particle at a site of the lattice Z, and - its 
absence. 

It is worthwhile giving an equivalent definition of the Markov measure without any 
connection to random processes (sec e.g. [25]). A measure fi on {0, 1}^ is Markov iff for 
any cylinder^ with the base where A and B two arbitrary finite binary sequences 

while 6 is a single binary letter, we have 

l,{[b])-i,{[AbC]) = i,{[Ab])-i,{[bC]). (3.1) 

Note that [25] claims that the Markov measures appear extremely rarely in cellular au- 
tomata systems and thus their existence for traffic maps turns out to be a surprise both 
in deterministic and random settings. 



3.1 Lattice traffic map 

Theorem 5 The deterministic lattice process tt{p = 1,v — l,r — 1/2, Z) possesses a 
Markov invariant measure ^, being a weighted sum of two measures of maximal entropy 
for shift-maps acting in opposite directions. Besides, there is a 1-parameter family of 
invariant measures {/xi,i,i/2,z,^}e such that E^^ ^ i/2Ze['^i] ~ 6- 

Note that the construction of Markov invariant measures which we shall elaborate for 
the analysis of the random setting in Section 4 cannot be applied here and we shall use 
a very different strategy related to the construction of measures of maximal entropy well 
known in theory of low-dimensional hyperbolic dynamical systems (see e.g. [10]). 

Let B := {0, 1}^ be the space of binary sequences and let B± be its two subspaces 

defined by the transition matrices M+ ^ i o ) '~ 1 1 ) ' ^^^'^^^ 

B+ := {6 e B : 6^ + 6^+1 ^ 2 Vi e Z}, 

B_ := {6 G B : 6^ + 6^+1 ^ Vi G Z}. 
Consider two shift-maps cr± : B ^ B defined as follows: 

{a+b)i := bi-i, {cr-b)i := 

i.e. these two maps shift a sequence by one position in opposite directions. Obviously 

(j±B± = B±. 



set of all sequences taking a given collection of values at a given collection of indices, called its base 



Denote by T : B — > B the map corresponding to the deterministic process tt{p — 
l,v = l,r = 1/2, Z). Recall that the density of a binary sequence 6 £ B is defined exactly 
as the density of a configuration of ones in b given in Section 1. 

Theorem 5 claims only the existence of Markov invariant measures but in fact we shall 
prove that the restriction of the dynamical system (T, B) to the Cantor sets B± possesses 
massive invariant measures in the sense that these measures are positive on each open 
subset. Before to give the proof of this result let us discuss reasons for the absence of 
massive invariant measures in the entire space. In Section 4 we shall show that in the 
true random case (0 < p < 1) for each a = poi G (0) 1) there is a nontrivial measure /xi^' 
(constructed in the proof of Theorem 7) whose value on the cylinder with the base [1100] 
is equal to 

/^aH[1100]) =PiPiiPioPoo -Pio)PioPoo- 
On the other hand, by (4.4) 

Pw = (1 -Poi)/(l -PPoi) ^ 1- 
Therefore for each poi we have 

/.«([1100])^0, 

and thus the limiting (as p — )■ 1) measure cannot be massive. This explains why in 
the deterministic case when p — 1 one should not expect the existence of true massive 
invariant measures. We conjecture that for each T-invariant probability measure r) there 
exists an open subset i? G B with rj{B) = 0. However at present we do not have a 
complete proof of this statement. 

Proof of Theorem 5. Let us start with a recipe of the construction of measures of 
maximal entropy for the shift-map cr : B^^ — )■ Bm with an irreducible transition matrix 
M = (mjj)ijg{o,i} (see e.g. [10]). For each n e denote by /x„ the probability measure 
uniformly distributed on points of period n of the map a. Then the measures /x„ weakly 
converge as n — > oo to a massive (on Bm) probability measure /la- The latter coincides 
with the unique invariant distribution for the stationary Markov chain with the transition 
probability matrix 

Pa := {mijmj/{\Mmi)), 

where \m is the maximal eigenvalue of the matrix M and (m^) is the corresponding 
eigenvector. The measure /Xg. represents the only invariant measure maximizing the metric 
entropy of the dynamical system (cr, B^) with the transition matrix M, which explains 
the reason for its name. Calculating the leading elements of the spectra of the matrices 
M± we get 

A+ = A_ =: A := (1 + V5)/2, m^ = m^ = l/X, mf = m^ = l-l/X. 

Despite the absence of a similar result for the dynamical system (T, B) we try to 
follow this recipe. A serious problem here is that we do not have good control over all 
time-periodic points of a given period for the map T. To overcome this difficulty we select 
a subset of time-periodic points with which we shall work. 

Observe that if a configuration 6 e B is spatially periodic then e Z+ T*b is spatially 
periodic with the same period. This fact follows immediately from the definition of the 



map T and is discussed in detail in [5]. Additionally in [5] it has been shown that if 
either (t"6 = h and h G B+ or a^h — b and b G B_ then this spatially periodic point is 
time-periodic for the map T with the same period n. Note that g{b) < 1/2 if 6 e B_|_ and 
g{b) > 1/2 if 6 e B_. 

The absence of non spatially periodic time-periodic points of the map T would com- 
plete the description of time-periodic points. Unfortunately this is not the case. The 
point is that the map is not one-to-one. Using this let us construct a sketch of a counter- 
example. Let b be spatially and time periodic for the map T with the same period n, i.e. 
bi = bi+n E Zi and T"6 = b. Thus a finite word 6-n+i, • • • , &-i) is a pre-image of the 
word (6o, ^1, • • • , bn-i) under the action of the map T". Since the map T is not bijective 
there exists another pre-image (&'_„, &'_n+i) ■ ■ ■ ■> ^-i) of the word {bg, bi, . . . , Simi- 
larly wc consider a pre-image of the word (6^, . . . , under the action of the map 
T" and denote it by {b'_2ni ^-2n+n • • • > ^-n-i)- Continuing this procedure we are getting a 
point &':=(..., b'_2i 60 &i, 62, • • •) ^ snch. that b' is no longer spatially periodic but 
= b'. 

Nevertheless we can use a spatially periodic part of time-periodic points for our con- 
struction. For each n e Z+ consider a probability measure /x* uniformly distributed on 
spatially and time periodic points of period n of the map T. As we already noted each 
spatially and time periodic point of the map T is spatially periodic point of one of the 
shift-maps (t±. Due to the symmetry of the motion of ones and zeros under the action 
of the map T the number C+ of n-periodic points of density less or equal to 1/2 differs 
at most by n from the number C~ of n-periodic points of density larger than 1/2, while 

arc of order r?/ (sec the calculation of the exponent A > below). Therefore one can 
represent the measure /x* as follows 

lA^ = (CnV: + c-/i;)/(c+ + C-), 

where 11^ are probability measures uniformly distributed on n-periodic points of the shift- 
maps (T±. Being uniformly distributed on time-periodic trajectories the measure /x* is 
T-invariant for each n. 

We already know that the measures //^ converge as n — > 00 to the corresponding 
measures of maximal entropy /x^ for the shift-maps (t±. Therefore 

lAx {l4 + /^a )/2 =: /^T- 

The weak massive property for the limit measure follows from the similar statement for 
the measures /x^ on the sets Bj.. Indeed by the construction for each cylinder on B_|_ its 
/xjr measure is positive, and the similar statement holds for each cylinder on B_ and the 
measure 

Note that according to our construction all three limit measures jix, fJ't T- invariant. 
Moreover, the equality of the leading eigenvalues A4. = A_ =: A := (1 + -\/5)/2 implies the 
coincidence of the corresponding metric entropies (being equal to In A) . 

It remains to construct for each g e [0, 1] the T-invariant measure //^ supported by 
the configurations of density g. This can be done using either Gibbsian reconstruction of 
the measures fi"^ or by an explicit representation in terms of Markov shifts. Let us discuss 
the latter approach. 

The action of a shift-map a is equivalent to the time-shift along realizations of the 
Markov chain with the transition probability matrix compatible with the corresponding 



binary transition matrix M. Recall that matrices A — (aij) and B — (bij) with nonnega- 
tive entries are com,patible if aijbij — implies aij + b^j = 0. 

A probability transition matrix compatible with the matrix M_|_ is written as P+ := 



bijection between the values of the parameter a and the set of particle densities g. Thus 
for each density g & [0, 1/2] we obtain a massive (spatially) Markov T-invariant measure 

Similarly one considers the matrix M_ which allows one to construct massive invariant 



Theorem 5 gives a complete recipe for the construction of nontrivial invariant measures 
for the process 7r(p — l,v — l,r — 1/2,Z). 

3.2 Existence of massive invariant measures in continuum 

In this Section we develop a special machinery extending the measures fi^ (constructed 
for the lattice case) first to nontrivial invariant measures of the deterministic process 
7r(p = l,v = l,r = 1/2, M) acting on the real line and then to deterministic processes 
7r(p = 1, V > 0, r > 0, M) with arbitrary v > 0, r > 0. 

Recall that in this more complicated 'continuous' setting a finite c?//m(ier with the base 
defined by a finite subset of integers / and a collection C := {Cjjjg/ of open intervals^ 
is the subset Cj^c '■= {x ^ X : Xi & Ci \/i & I}. We endow the space of admissible 
configurations X with the cr-algebra B generated by the finite cylinders defining a topology 
in this space. 

By Lemma 1(1) the measures /ig — fig^r are 7r(p — l,v — l,r — 1/2, R)-invariant. 
For each r' > applying the affine transformation (2.1), obtained in Lemma 1(3), to the 
measure we arc getting a new probability measure fi^'^r' with g' := g/{l + 2{l/2 — r')g), 
which is tt{p = 1,v = 1, r', M)-invariant. Since g takes all values from the interval [0, 1] 
the new variable g' takes all values from the interval [0, l/(2r')]. 

Making yet another spatial change of variables z — >■ vz + w with parameters v > 
0, w >0 we are obtaining from the measure n^'^r' a two-parameter family of probability 
measures ^{p = 1, g", r", v, w) supported on configurations of balls of radius r" := vr' and 
having density g" := g'/v. Applying again Lemma 1(4) we see that for each g",r',v,w 
the measure /x(p = 1, g", r", v, w) is 7i{p = 1, v, r", M)-invariant. 

On the other hand, these measures cannot be massive, since they are supported by 
very thin sets. To construct a massive (albeit non-crgodic) invariant measure from this 
family we consider a measure n{p = 1, g", r", v) having marginals ii{p = 1, g", r", v, w) on 
sub-lattices fZ + w and uniformly distributed with respect to the parameter w G [0,v). 

This finishes the proof of the existence of massive invariant measures for the deter- 
ministic traffic map in continuum. 

The construction of the massive invariant measure in the lattice case (R := Z) with 
general v G Z+ is similar except that w G {0, 1, . . . ,v — 1} and r" = 1/2 or r" = to be 
able to work with point particles in lattice setting. 




with a single parameter < a < 1. The function g :— a/{l + a) defines a 



measures /ig with ^ G (1/2, 1]. 



□ 



^In general the cylinder Ci^c might be empty for nonempty sets I, C. 



3.3 Entropy 



To finalize the description of ergodic properties of the processes under study in the pure 
deterministic setting we show (following the approach developed in [7]) that these pro- 
cesses are strongly chaotic. Our dynamical system is defined by a deterministic map 
T^, : X — )■ X from the set of admissible configurations into itself. Our aim is to show that 
the topological entropy of this map is infinite.^ 

We refer the reader to [26, 27] for detailed definitions of the topological and metric 
entropies for deterministic dynamical systems and their properties that we use here. To 
avoid difficulties related to the non-compactness of the phase space we define the topo- 
logical entropy of a map Ty (notation /itop(^i))) as the supremum of metric entropies of 
this map taken over all probabilistic invariant measures (compare to the conventional 
definition of the topological entropy and its properties discussed, e.g. in [27]). 

We start the analysis with the action of a shift-map in continuum : X ^ X defined 

as 

(cr^x)j :— Xi-\- V i & Z, x & X, v > 0. 
Lemma 3 The topological entropy of the shift-map in continuum is infinite. 

Proof. The continuity of the shift-map in continuum in the topology induced by the 
(T-algebra B generated by finite cylinders is implied by the fact that a preimage of a finite 
cylinder under the action of is again a finite cylinder. 

The idea of the calculation of the entropy is to construct an invariant subset of X on 
which the map is isomorphic to the full shift-map in the space of sequences with a 
countable alphabet. The result follows from the observation that the topological entropy 
of the full shift-map cr^"^^ with the alphabet consisting of n elements is equal to In n (see, 
e.g. [26, 27]). 

Let a := {ai}i(zz_^_ with ctj G {0,v) and let a*^"^ := {ai}^^^. Consider a sequence of 
subsets X(") C X consisting of all configurations x E X satisfying the condition 

X2k e vZ, X2k+i e + oi^"^ VA; e Z. 

Then X^") is (7„-invariant and the restriction av\X^'^^ is isomorphic to the full shift- map 
cr^") with the alphabet consisting of n elements {cj} of type aj := {[0, ctj), [oj, v)}, i.e. 
each element is represented by a pair of neighboring intervals. Therefore the topological 
entropy of a^"'^ is equal to Inn ^^—^ oo. □ 

Theorem 6 The topological entropy of the traffic map in continuum M is infinite. 

Proof. The traffic map is continuous in the topology induced by the o"-algebra B 
generated by finite cylinders by the same argument as in the case of the shift-map. 
Observe that the subset 

X>, {xeX : Ai{x) >v Vi G Z} 

■^Normally one says that a map is chaotic if its topological entropy is positive, so infinite value of the 
entropy indicates a very high level of chaoticity. 



of the set of admissible configurations is T^-invariant. Therefore 

and for our purposes it is enough to show that the latter is infinite. On the other hand, 
by the definition of the map we have Tt,|X>,u = o',u|X>^. 

One cannot apply the result of Lemma 3 directly because in the case under consider- 
ation the gaps between particles are greater or equal to v by the construction, while in 
the proof of Lemma 3 the gaps did not exceed v. To this end one sets ctj e {v,2v) and 
modifies the definition of X^") as follows: 

^^"^ := {X2k e 3vZ, X2k+i e X2k + a^"^ VA; e Z}. 

Consider the the alphabet A^"'^ with elements of type aj := {[0,Q;j), [ai,3v)}, ai E a^'^\ 
Then the 3-d power of the map T^|X>^ is isomorphic to the full shift-map cr*^"^ with the 
alphabet ^4^"^^ Using that 

3/itop(7;|^>.) = htop{{T,\X>,f) = /itop((J^"^) = Inn 
we get the result. □ 



4 Ergodic properties: stochastic case (0 < p < 1) 

4.1 The simplest lattice process 7r{p, v = l,r = 1/2, Z) 

We start with the construction of a nontrivial Markov invariant measure for the simplest 
lattice TASEP. 

Lemma 4 Let P ^ ^ be a probability matrix with positive entries and the left 

leading normalized eigenvector {po,Pi). Then the measure /i on X, defined on cylinders 
by the relation 

n-l 

/i([ai,a2, . . . ,a„]) := Pa^ JJ^ e{0,i}, (4.1) 

i=l 

is invariant with respect to the process Tr{p, v = l,r — 1/2,Z) iff 

PooPii = (1 - p)PioPoi- (4.2) 

Proof. For a Markov chain n acting on a space of binary sequences we say that a 
cylinder S is a (partial) pre-image of a cylinder A if the probability Pr(7r(S) — A) > 0. 
We need to check that for any cylinder its measure is equal to the sum of measures of all 
its pre-images multiplied by the corresponding transition probabilities. The proof follows 
by induction on the cylinder's length n. 

We start with cylinders of length 1 and 2. The following tables show all pre-images of 
cylinders of length 1 and 2. The left column corresponds to cylinders and their pre-images, 
while the right column shows stationary probabilities for the cylinder (the 1st line) and 



stationary probabilities for its pre-images multiplied by the transition probabilities. Here 
p :— 1 — p is the probability that the jump does not take place. 
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PiPn 
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PPiPioPoo 
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PPPiPwPoi 


1010 


P^pApoi 


1011 


PPiPioPoiPii 


01 


PoPoi = PlPlO 




10 


PPiPw 






0010 


PPoPooPoiPio 


= PPiPwPoo 


0011 


PoPooPoiPii = 


PPipIoPoi 


1010 


P^PiPloPm = 


PPiPwPooPu 


1011 


PPiPioPoiPii 







In the first two cases the equivalence of the probability in the 1st line to the sum of 
other probabilities is trivial. In the 3d case assuming that ppiQ 7^ we get the following 
condition for the equivalence 

1 =p + ppn + ppoo + P^PioPoi =^ 1 = Poo + Pii + PPioPoi, (4.3) 
while in the 4th case assuming that piPn 7^ we have 

1 =PPio +P11 +PPP?oPoi/Pii +PP10P01 =^ PooPii =PPioPoi- 

The last calculation proves that the assumption (4.2) is necessary. A direct calculation 
shows that despite appearances the relation (4.3) is equivalent to (4.2). 

The checking of the cases 5 and 6 can be done similarly using the property poPoi = PiPio 
(which we already used in the corresponding tables). 

Assume now that the claim is already proven for all cylinders with bases less or equal 
to n > 1. To reduce the case of length n + 1 to n observe that the only difference between 
smaller cylinder's length to the larger one consists of an additional letter or 1 at the 
right end of the cylinder's base. Therefore it is enough to show that the corresponding 
probability changes by Pab, where a is the last letter in the shorter cylinder and b is the 
new letter. To prove this we consider all 4 possibihties: 

1 ^ 10,1 ^ 11,0 ^ 00,0 ^ 01. 

The 2nd and 3d situations are relatively simple: 

1 ^ 11 : 10 ^ 110, 11 ^ 111 and 10 ^ 1010 + 1011 

PPiPio(pPioPoi +P01P11) = PPiPio(pooPii +P01P11) =PPiPioPii- 

^ 00 : 00 ^ 000, 10 ^ 100 and .10 ^ 0010 + 1010 
(here and in the sequel "." stands for an arbitrary symbol) 



PPiPio(pioPoo +PP01P10) = PPiPio(pioPoo +P00P11) = PPiPioPoo- 



Note that we use that the transition probabilities depend only on the previous letter. 



In the 1st situation we proceed as follows 

1^10: 10 + 11^10 + 110, 10^100 + 1010. 

Denoting the product of the measure of a cylinder A and the corresponding transition 
probability by K{A) we get 

K{[10]) + K{[11]) =PiipPio+Pn) = piPio - PPio + Pn) = Pi{^ - PPio) 
K{[10]) + K{[im]) = pipioip + ppu) = piPio(l -P + PPii) 

= PiPio(l - P(l - Pn)) = PiPio(l - PPio)- 

Thus the 2nd sum differs from the 1st by the desired multiplier piQ. 

The 4th situation is a bit more complicated. In this case we split the 3d pre-image 
.10 of into two parts 010 and 110 giving contributions 

X([010]) = ppoPoiPio, i^([110]) = PPiPiiPio 

respectively (observe that ppoPoiPio + PPiPiiPio — PPiPio), ^-nd the 1st pre-image 10 of 01 
also splits into two parts 010 and 110 giving contributions 

i^([010]) = ppoPoiPio = PPipIo: ^([110]) = PPiPnPio 

respectively. Then we gather them as follows: 

00 + 010 ^ 010 + 0010 + 0011, 10 + 110 ^ 110 + 1010 + 1011. 

The following simple calculation checks the correctness of this construction: 

K{[00]) + K{[010]) = poipoo + PPoiPio) = Po(l - Pii) = PoPio 
K([0010]) + i^([0011]) + K([010]) = PipUpPoo + PPoi + p) 

= PiPio = PoPoiPio = (PoPio)Poi- 

□ 



Theorem 7 The process 7r{p,v — l,r — 1/2, Z) Vv e Z\,p,g e (0,1) possesses a 1- 
parameter family of probabilistic invariant measures {Afp,i/2,i,z,e}e positive on each open 
set and supported by the set of configurations of density g. 

Proof. Observe that the probability measure /j, constructed in Lemma 4 due to the 
additional condition (4.2) can be parametrized by a single parameter a :— poi ^ (0, 1). 
Prom (4.2) we obtain 

Pio = (l-a)/(l-pa). (4.4) 

In view of the equality between the density of particles g and the stationary probability 
of ones pi under the Markov shift with the transition matrix P we have 



g = a/{a + pio) = a(l — pa)/ (1 — pa^). 



Solving the last equality with respect to the parameter a, we get 

^ 1 - yi - 4pg(l - g) 
2p{l - g) 



(4.5) 



Therefore, for a given p the constructed family of measures {/^^^l^ is uniquely indexed 

by the density g. 

The property that the measures {/U^^''}^ are massive follows immediately from the 
observation that the measure of an arbitrary finite cylinder is positive. □ 



4.2 Existence of massive invariant measures in continuum 

Existence of massive 7r(0<p<l,f = l,r = 1/2, Z) -invariant measures iJ,f^ for each g is 
already proven in Section 4.1. 

By Lemma 1(1) the measures /j,^^ are 7r{p,v — l,r — 1/2, R)-invariant. For each 
r' > applying the affine transformation (2.1), obtained in Lemma 1(3), to the measure 
//^'* we are getting a new probability measure jig/y with g' :— gj (1 + 2(1/2 — r')^), which 
is 7r{p,v = 1, r', M)-invariant. Since g takes all values from the interval [0,1] the new 
variable g' takes all values from the interval [0, l/(2r')]. 

Making yet another spatial change of variables z ^ vz + w with parameters v > 
0, ty > we are obtaining from the measure /ig/y a two-parameter family of probability 
measures iJ,{p, g",r",v,w) supported on configurations of balls of radius r" := vr' and 
having density g" := g'/v. Applying again Lemma 1(4) we see that for each g",r',v,w 
the measure /x(p, g", r", v, w) is 7i{p, v, r", M)-invariant. 

On the other hand, these measures cannot be massive, since they are supported by 
very thin sets. To construct a massive (albeit non-ergodic) invariant measure from this 
family we consider a measure fi{p, g" ,r" ,v) having marginals ^{p, g" ,r" ,v,w) on sub- 
lattices + w and uniformly distributed with respect to the parameter w G [0,f). 

The construction of the massive invariant measure for 7r{p,v, r, Z) with general v € 
is exactly the same except that w e {0, 1, . . . , t; — 1} and r" = 1/2 or r" = 0. 



4.3 Stochastic stability 

Despite very substantial differences in the constructions of nontrivial invariant measures in 
deterministic and random situations the structure of the invariant measure in both cases 
is described in terms of 2 x 2 Markov matrices {p^^}- Therefore to prove the stochastic 
stabihty one only needs to check that 

To this end one uses exphcit formulas for the entries of the Markov matrices. Namely by 
(4.4) 

pSS = (l-a)/(l-pa)^l, 

which coincides with the corresponding entry in the deterministic case, while p^f^ 0. 
Passing to the limit as p — > 1 one gets the corresponding limit relations for two other 



entries pocPoi, which define uniquely the Markov measure in the lattice deterministic 
setting with v = 1. Thus the nontrivial invariant measures for the process 7r(p — l,v — 
l,r = 1/2, Z) are stochastically stable. 

Recall now that the nontrivial invariant measures for the exclusion type processes 
in continuum 7r(p,t', r,M) were constructed in Sections 3.2,4.2 in three steps through 
changes of variables which do not depend on the choice of the variable p. Therefore this 
construction withstand the limit transition as p ^ 1, which proves the stochastic stability 
for the general process n{p = l,f,r, M). Similarly one proves the stochastic stability for 
the lattice processes with long jumps 7r(p = 1, v e Z+, r e {0, 1/2}, Z). 

Let me note that the Markov measures are unique among stochastically stable trans- 
lationally invariant measures, but there are non-translationally invariant ones. To demon- 
strate this consider a 5-measure fi supported by a single configuration x :—{... ,0,Xo — 
0,Xi = 1, 1, . . .). Then this configuration is a fixed point for the process 7r{p,v — l,r — 
1/2, Ia) for each < p < 1 and hence the measure /i is stochastically stable. 



5 Average velocities 

Explicit expression for the invariant measure n allows one to derive by Birkhoff 's ergodic 
theorem the formula for the average velocity for /i-a.a. initial configurations. In what 
follows we are interested in a much more broad setting of all particle configurations having 
well defined densities. This extension can be justified as follows. Theorem 3 shows that 
for all particle configurations of a given density the average velocity is the same (or for 
all of these configurations the average velocity is not well defined) . Therefore it is enough 
to calculate this statistics for a most suitable single initial configuration of given density. 
Choosing a configuration typical with respect to the invariant measure fi we achieve this 
goal. Since this construction does not depend on specific properties of the invariant 
measure we shall not repeat this argument in further calculations. 



5.1 Lattice TASEP with v = l 

Let us use the constructed Markov measures to calculate the average velocities for the 
process n{p,v,r = 1/2, Z). 

Lemma 5 V{g) = (1 - ^1 - 4p^(l - g))/{2g) for 7r{p <l,v,r = 1/2, Z). 

Proof. For a given configuration of particles, a particle may move iff the next site of 
the lattice is not occupied, i.e. only in the situation 10. Therefore the average velocity 
is equal to the stationary probability of the jump to the right, which in turn is equal to 
ppiQ. The representation of the Markov invariant measure obtained earlier immediately 
gives the formula for the average velocity 

V{0,p,l,l/2)=ppio=pil-a)/{l-pa) E [0,p]. 

Substituting the value a = a{g) according to the formula (4.5) we get 

V{g,p, v^l,r^ 1/2) = [l - ^1 - 4p^(l - g)] /{2g). (5.1) 

□ 



5.2 Process in continuum 



By means of results of Theorem 3 and Lemma 1 the relation (5.1) can be transfered to 
the processes of the 3d type 7r(p, T;,r, R) without changes. Note that the formula (5.1) 
was already known in physics publications for the lattice processes of type 1 with v — 1 
(see [22]). 

It is important to say that the naive transition from v = 1 to v > 1 directly in the 
class of lattice processes (using the invariance of sub-lattices with step multiple to v) is 
impossible (or rather so we can study only low-density < 1/v configurations). Instead, 
we use the self-similarity of processes of type 3 acting in a continuous space. 

Lemma 6 Let the n{p,v,r,M) type processes a;*,i* with parameters r = r = 0, v > 
0, V — uv, p = p E (0, 1) having initial configurations — ux^ with some u > be 
statically coupled. Assume also that g{x) and V{x) be well defined. Then 

g{x) = g{x)/u, V{x) = uV{x). (5.2) 
Proof. By a straightforward calculation. □ 

Applying these similarity transformations to the special case described in the relation 
(5.1) we obtain the general formula for the average velocity (1.3). 

CoroIIeiry 7 

V{g,p, V > 0,r — 0) — vV{vg,p, v = l,r = 1/2) 



1 + vg — -\/(T+^^^P"-^4pv^ 
mm{v, 1/g). 

Since the last term above coincides with the average velocity for the deterministic 
traffic map (see [7] ) , the limit transition in the last relation shows that the average velocity 

is stochastically stable. 

Indeed, under a spatial change of variables 1 — )■ f and the transition from the config- 
uration of balls of radius r = 1/2 to the configuration of point-particles (i.e. r = 0) with 
the same sequence of gaps we get 

Q^q/C^ - 2r^) = g/{l - g), 

hence g = g/{l + g) and 



V{g,p,v = l,r = 0) = ^ — 



1 - ipj^ l - ^ 

l + g-\/(l + g)^-4pg 
2g 



Therefore by Lemma 6 



V{g,p,v > 0,r = 0) 



= vV{vg.p. V = l.r = 1/2) 



1 1 + vg — ^/(l + vg)"^ — ApvQ 




l + VQ- 1/(1 + vqY - Apvg p^i 



mm{v, l/g). 



2g 



This finishes the proof of Theorem 2. 

5.3 Lattice exclusion processes with long jumps 

The results for 7f{p, v, r, M) with arbitrary v G Z+ are transferred directly by Lemma 1 
back to the lattice cases. Indeed, Theorem 3 shows that it is enough to derive the formula 
for the average velocity for a specially chosen initial configuration of a given density. On 
the other hand, by Lemma 1(1 and 2) a reahzation of the processes 7r{p,v,r — 1/2, Z) 
and 7r(p, v, r = 0, Z) starting from certain configurations coincide with a realization of the 
7i{p,v,r,M) process statically coupled to the lattice process and starting from the same 
initial configuration. 

It is worth noting that a naive application of the property described in Lemma 1(4) 
seems to extend the results about the process 7r(p, v = 1, r = 1/2, Z) directly to n{p, v > 
l,r — 1/2, Z) restricting the latter process to invariant sub- lattices vZ+w, w — 0,1, . . . ,v— 
1. A close look shows that this is indeed the case but only for configurations of low den- 
sity g < 1/v, since otherwise particles from the same configuration located at different 
sub- lattices will interact. 

5.4 Heterogeneous particles (of different sizes) 

Thinking about the processes under consideration as models of traffic ffows it is reasonable 
to take into account that vehicles need not to be of the same size. From this point of view 
we consider an exclusion type process in continuum with particle configurations consisting 
of balls with varying sizes, i.e. the radius of the i-th ball is equal to > 0. Our aim is 
to show that the dependence of the average velocity on density in this case can be easily 
obtained from the corresponding result for the case of balls of the same radius. 

Let a; be a bi-infinite admissible configuration of particles represented by balls of 
radiuses n > with the average value 



and centered at points Xi e M. The notion of admissibility and the law of dynamics should 
be slightly rewritten (in comparison to the homogeneous case): 




n-l 



r := lim — 

n—i-oo n 



i=0 



Xi + Ti Xi^i — rj_|_i. 




'i+l 



Tj+i} with probability p 

with probability 1 — p 



Theorem 8 Let the process defined above and the Tr{p,v,'r,'R) process be statically 
coupled and let r — r. Then the average velocities of these processes coincide. 

Proof. A simple generalization of the affine conjugation between configurations of 
different ball's sizes introduced in Lemma 1 allows one to make the bijection between the 
configurations x*, f*. Indeed, consider an affine map defined by the relation 

i 

{(p{x))i :^Xi-2j2irj-r). 

j=0 



Then 



ft ft 

n-1 



- 4) + I E(^^- - - ^(^0 - r) 
l/g{x^) + 2{r-r). 



n n ^ — ' n 

3=0 



Thus ^(x*) = ^(x *) if r = r. 

Now we are ready to calculate the average velocity V{x): 

V{x,i,t) = \{x\- x^) = \{{^-\x'))^ - {ip{x%) 



r ' t 

j=0 



t—>-oo 



V{x) + 2{f-f), 

which proves our claim. □ 



5.5 Heterogeneous space 

So far we have considered only exclusion processes acting on homogeneous spaces. In [8] 
we introduced and studied a modification of the deterministic version of the exclusion 
process in continuum which takes into account the presence of static obstacles (traffic 
lights) for the motion of point particles (i.e. r — 0). Fix an arbitrary point-particle 
configuration z = {zj)j^z G X(0, M) whose elements correspond to positions of obstacles. 
Then the formula (1.1) can be rewritten as follows: 

^t+i ^ rmin{x* + V, Zj^^^t^} with probability p 
* \x\ with probabihty 1 — p ' 

where j{xf) :— min{/c e Z : x\ < Zk}. Thus the "obstacles" suspend the movement of 
particles, taking into account the time necessary to overtake an obstacle. 

For a given f > and a configuration of obstacles z denote by z the extended con- 
figuration of obstacles obtained by inserting between each pair of entries Zi.ZiJ^i new 
\_{zi^i — Zi)/v\ 'virtual' obstacles at distances v between them starting from the point Zi. 
Here \u\ stands for the integer part of the number u. 



Theorem 9 For given v > Q < p < 1 and any configurations x,z E X for which the 
densities g{x), g{z) are well defined 



Note that the average velocity in the above formula does not depend explicitly on 
the local velocity v, however the latter is included to the construction of the extended 
configuration z, in particular g{z) > l/v. It is interesting to note also that in [8] it was 
shown that under more general setting with the non-degenerate distribution of random 
iid local velocities vj the average particle velocity may not exist. 

The difficulty of the analysis here is is that the inhomogeneity of the space in which 
the collective random walk takes place (the presence of obstacles), generally does not 
permit the existence of invariant measures.^ Therefore, the main step of the approach we 
used - the construction of a massive invariant measure is impossible. The complete proof 
of this result needs a modification of the dynamical coupling construction elaborated in 
[7, 8] and will be discussed elsewhere. Here we only give an idea of the proof. 

Using the technique developed in [8] for the deterministic version of the problem it 
can be shown that the calculation of average velocities V{x, z, v,p) can be reduced to the 
analysis of a Markov process of type 2, acting (in contrast to the already-studied setting) 
on a inhomogeneous lattice R := z. The existence of the density of the configuration z 
allows one to transfer the results obtained for the conventional integer lattice Z to the 
inhomogeneous case under consideration. This completes the construction. 



^For the existence of invariant measures one needs at least the condition of stationarity for the con- 
figurations of obstacles z. 



V{x,z) 



g{x) + g{z) - ^/igix) + g{z)y - 4:pg{x)g{z) 




(5.4) 
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